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Abstract. We study the dynamics of stellar wind bubbles around hydrogen-deficient stars using numerical simulations with 
time- and ion dependent cooling. We consider two types of hydrogen-deficient stars, massive WR stars, producing Ring 
Nebulae, and low mass [WR] stars, producing Planetary Nebulae. We show that for the Planetary Nebulae, the different cooling 
properties of the hydrogen-deficient wind lead to a later transition from momentum- to energy-driven flow, which could explain 
the observed turbulence of these nebulae. We find that Ring Nebulae should all be energy-driven, and show how comparing 
the bubble's momentum and kinetic energy to the input wind momentum and kinetic energy, can give misleading information 
about the dynamics of the bubble. 
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1. Introduction 

Both high and low mass stars can under certain circum- 
stances reduce the hydrogen content of their atmospheres. In 
most cases this leads to the so-called Wolf-Rayet (WR) phe- 
nomenon, i.e. a dense fast wind starting below the photosphere, 
which produce s a Wolf-Rayet spectr um, dominated by bright 
emission lines ( A.bbott & Conti 1987 ). 

Traditionally most attention was given to the massive WR 
stars, but the last ten years their lower mass c ousins, the [WR] 
stars have been studied in more detail (see e.g. Koesterke 2001 , 
for a review). 

The winds from [WR] and WR stars produce stellar wind 
bubbles (SWBs) in their environment. In the case of WR stars 
they are called Ring Nebulae (RNe), in the case of [WR] stars. 
Planetary Nebulae (PNe). PNe also form around H-rich cen- 
tral stars, so [WR] stars constitute a subgroup among central 
stars of PNe. Approximately 7% of central st ars is estimated to 
be [WR], the rest being H-rich ( Gorny 2001 ) (with the excep- 
tion of so-called weak emission line stars or wels which appear 
to be H-poor without showing the WR phenomenon). All cen- 
tral stars are considered to be in the same evolutionary phase, 
namely the post-AGB phase, where the [WR] have changed 
their atmospheric abundances through a timely thermal pulse 
( Herwi j 2001 ). The existence of two different groups of cen- 
tral stars suggests that a comparison between the two could be 
interesting. 

In contrast, WR stars are thought to be a phase in the evolu- 
tion of most stars with Zero Age Main Sequence masses higher 
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than ^ 25 Mq, and there is no class of H-rich stars in an equiv- 
alent evolutionary stage. Massive stars lose mass already on 
the Main Sequence, followed by a slower wind when the star 
moves over to the red part of the Hertzsprung-RusseU diagram, 
and finally leading to the WR wind. This leads to a whole se- 
ries of interactions between the wind phases (see e.g. Garcia- 
Segura et al. 1996b). The more complicated environments and 
probably also the dumpiness of the actual winds make that the 
RNe are mostly irregular and filamentary, lacking the overall 
symmetries found in PNe. 

In this paper we investigate whether the fact that the winds 
from WR and [WR] stars are H-poor changes the dynamics 
of their SWBs. In the case of the [WR] stars this is relevant 
because we can compare the PNe between the H-rich and H- 
poor central stars. In the case of the WR stars this is relevant 
because it has been suggested that even at high wind velocities 
their SWBs can be strongly cooling. 

The layout of the paper is as follows. The effects of cool- 
ing on SWBs are outlined in Sect. 2. We investigate the effects 
of WR winds using numerical hydrodynamic models with de- 
tailed cooling, described in Sect. 3. Section 4 and 5 contain the 
results of the simulations for PNe and RNe, respectively. We 
discuss these results further in Sect. 6 and sum up the conclu- 
sions in Sect. 7. 

2. Stellar wind bubbles 

Stellar wind bubbles (SWBs) expanding into a dense environ- 
ment (as is the case in both RNe and PNe), come in two basics 
types: radiative (also called momentum-driven or momentum- 
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conserving) and non-radiative (also called adiabatic or energy- 
driven/energy-conserving). Whether the SWB is radiative or 
non-radiative depends on the cooling at the wind shock (or in- 
ner shock). This is a well-known phenomenon, see for example 
Dyson & Wiliiams| ( |l997| ) or |Koo & McKe"^ ( [l992^ ,|b|). 

Since cooling strongly depends on the metallicity, it is ex- 
pected that the radiative cooling rate of the metal-rich WR 
winds will be orders of magnitude above normal ('solar abun- 
dance') cooling rates. This could then lead to SWBs remaining 
radiative when they would normally be non-radiative. 

Whether a SWB is radiative or non-radiative depends on 
three time scales (Koo & McKee 1992b| ): the crossing time for 
the free wind Across = ^sw/^'w, the age of the bubble t, and the 
cooling time tcooi- The cooling time can be found from 



^cool 



Ppre 



(1) 



where Ci is a constant which describes the cooling, which is 
normally taken to be 6.0 x 10^"^^ g cm^^ s'*. If tcooi ^ icross 
the bubble is radiative, if t <C icooi the bubble is non-radiative, 
and if tcross ^ icooi «C t cooling does affect the shocked wind, 
but the bubble is still mostly filled with hot shocked gas, and 
the bubble is partially radiative. 

The application of these expressions lead to evolutionary 
sequenc es for stellar wind bubbles, as explored by Koo & 
McKee (19921:) for the most general case, by Garcia-Segura 



& Mac Low (1995) for RNe, and by Kahn & Breitschwerdt 



(1990) for PNe. Here we limit ourselves to the interaction of 



constant or evolving winds interacting with a previous slow 
wind phase, the generic model for both PNe and RNe. When 
a volume of hot shocked wind material is present, we will call 
the bubble energy-driven, if not, we will call it momentum- 
driven. This means that a partially radiative bubble will also be 
considered to be energy-driven. 

3. Numerical method 

We study the formation of SWBs using a numerical hydro- 
dynamics code coupled to a detailed atomic physics mod- 
ule, which calculates the ionization/recombination and cool- 
ing for the ions included. The hydrodynamic Euler equations 
are solved in one spatial dimension using the Roe-solver (Roa 



1981; Mellema et al. 1991). The spatial coordinate is radius. 



The Roe solver solves the equations conservatively, and with 
second order accuracy in both space and time. The ionic/atomic 
concentrations are passively ad vected using the approach from 
CLAWPACK dLeVequel ^002| ). This advect ion scheme was 
tested using the test problems described in ^lewa & Miiller 



( |1999[ ) and found to perform as well as the approach described 
by those authors. 

3. 1. Atomic physics module 

Radiative cooling is due to excitation and ionization of 
molecules, atoms and ions, and therefore depends on the com- 
position of the gas. Most numerical simulations which include 
cooling, use a so-called cooling curve. This is a curve which 
gives cooling as a function of temperature. The first cooling 



curves were constructed in the sixties, and the one which be- 
came standard were done around 1970 by Cox and collabora- 
tors ( |Cox & Tucker! |1969| ; |Cox & Daltabuil| |1971| ; Raymond 
et al. 1976). These curves are calculated assuming collisional 



ionization equilibrium to fix the ion concentrations. Although 
the original publications contained cooling curves for individ- 
ual elements (but not ions), the most common use of them has 
been in the form of one total cooling curve for standard abun- 
dances, such as the one from Dalgarno & McCray (1972). For 
cosmological applications, H and He only cooling curves have 
been constructed; Sutherland & Dopita (1993) present cooling 



curves for a range of metallicities. All of the above, as well as 
our method below, assume an optically thin plasma, in which 
all of the radiation escapes. 

Primarily motivated by astrophysical situations in which 
photo-ionization (and not collisional ionization) determines the 
ion concentrations, we have developed numerical approaches 
in which cooling is calculated for each ion separately and time- 
dependently, using the local tem perature and densi t y. The first 
version of th i s was described by Frank & Mellema ( 1994) and 
Balick et al. { 1993 ), and a newer version by |Raga et al~ ( 1997 ) 
and iMellema et al. (|1998|). 



Here we use an improved version of the same approach as 



in Raga et al. (1997), and VIellema et al. (1998). The improve- 
ments are that we now allow for position-dependent elemental 
abundances, which is crucial for the problem studied here. We 
also now include all ions of the following elements: H, He, C, 
N, O and Ne. Collisional excitation for the highest ionization 
stages were taken from Gaetz & Salpeter (1983). We have also 
included free-free and free-bound cooling to allow an extension 
to 10® K using the method of Gronenschild & Mewe (1978). 
Free-free Gaunt factors are from ^arzas & Latter ^ 961), and 
relativistic free-free corrections from poilld U980) 



Furthermore, we dropped the special treatment of H and 
He, which are now treated in the same way as the other el- 
ements. Since abundances can vary, we let all elements con- 
tribute to the electron density. Cooling due to collisional ion- 
ization is currently not included. This will underestimate the 
cooling at lower shock velocities. The module as used for the 
results in this paper does not contain photo-ionization, since we 
are primarily interested in the effects of the H-deficient winds. 

Our module does not include any heavier elements than Ne. 
At the high temperatures where the elements included become 
completely ionized, heavier elements can start to contribute to 
the cooling. We therefore underestimate the radiative cooling in 
the temperature regime logig T = 6 to 6.8, beyond which the 



free-free starts to dominate (see e.g. [Gaetz & Salpetei 1983) 
Whether or not this is a serious underestimate of the cooling 
depends on the relative abundances of the cooling elements. In 
our simulations it will mostly affect the solar abundances case. 
We will comment further on this below. 

In Fig. |l] we show three equilibrium cooling curves calcu- 
lated with the DORIC atomic physics module. One for solar 
abundances, one for primordial abundances, and one for [WC] 
abundances. These curves illustrate the dramatic effect of abun- 
dances on the cooling. Note that the high cooling rate below 
lO'' K is entirely due to CII, since this ion is assumed never to 
recombine (due to a CI ionizing background). 
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Equilibrium cooling curves 



Temperature (K) 

Fig. 1. Collisional equilibrium cooling curves calculated with 
the DORIC module. Solid line: solar abundances, dashed line: 
metal free, dotted line: [WC] abundances. The curves are not 
used in any of the calculations, but are shown to illustrate the 
effect of metallicity and the capability of DORIC. 

In the calculations below, we do not use these curves, but 
instead calculate the concentrations of ions (and therefore the 
radiative cooling) time-dependently. Especially in low density 
areas the ionic concentrations can deviate from the equilibrium 
values. We comment on this in Sect. 5. 

4. Models for PNe 

The properties of [WR] stars have been well studied, and the re- 
sults from two inde pendent groups ag ree well. For a n over view 
see the reviews of Koesterke ( 2001 ) and paman ^ ( |2002| ). All 
[WR] stars are carbon-rich, so from now on we refer to them 
as [WC] stars. What appears to be an evolutionary sequence of 
[WC]late, to [WC]eai-ly, to the PG 1 159 stars, is well sampled, 
and the abundances are more or less the same all along this 
sequence; a previous discrepancy has now been resolved, (De 
Marco, private communication). 

The stellar atmosphere analysis yields mass loss rates and 
velocities. The results show that the wind velocities are similar 
to the ones from normal central stars, but the mass loss rates 
are ~10 times higher. Therefore, we will use a higher mass 
loss rate in our [WC] model. 

To study the effects of H-deficient [WC]-type winds on 
the formation of PNe, we use the approach from Kahn & 
Breitschwerdt ( [1990[ ) and pwarkadas & BaHck| ( |1998| ). The 
evolution of the stellar wind in the post- AGE phase is not well 
known, not for normal, nor for [WC] central stars. Hence a sim- 
ple power law behaviour is postulated. Following Dwarkadas & 



Balick (1998) the fast wind evolves according to 



Wfw = Wfw(O) ( 1 + " 



(2) 



The initial velocity Wfw(O) is taken to be 25 km s~^, a = 
1.5, and t is chosen such that at < = 3000 years, Vf^ = 
2000 km s^^. The mass loss changes so that the product 
Aifwt^fw remains constant. This is expected for a radiation- 
driven wind (of constant efficiency). 



The environment is taken be a slow AGB wind character- 
ized by a constant mass loss rate of 2 x 10^^ Mq yr^^ and a 
velocity of 10 km s^^. Its abundances are solar. 

We calculated the evolution of this system for two cases, 
one for a fast wind with normal, solar abundances, and one for 
a fast wind with [WC] abundances. These abundances are as 
listed in Table 1 . The wind speed evolution is the same in both 
models, but the mass loss rate is taken to be 10 times higher in 
the [WC] case. The initial values (at a velocity of 25 km s~^) 
are 10"^ M© yr'^ for the [WC] star, and 10^^ Mq yr^^ for 
the H-rich star. As pointed out above, the mass loss rate drops 
as the velocity increases. At the end of the calculation when 
the wind has a velocity of 2000 km s^^, the mass loss rates are 
10^^ and 10~^ Mq yr~\ respectively. 

Figure || shows the evolution of the logarithm of the den- 
sity as function of radius and time. Light shades correspond to 
higher densities. In this representation the evolution of the shell 
is easily traced. In both cases the bubble is initially momentum- 
driven and the inner shock and contact discontinuity lie close to 
each other, at the inner edge of the swept up shell. In the case of 
a fast wind with solar abundances, the inner shock detaches it- 
self from the contact discontinuity at a time of 450 years, when 
the stellar wind has a velocity of 150 km s^^, and the shell has 
reached a radius of 7.5 x 10"^ pc. This is in line with what 
was found by previous authors ( ^<^ahn & Breitschwerdt 1990 



pwarkadas & Balick) |1998| ). In the case of the [WC] star, the 



transition happens at a much later time, 1100 years after the 
AGB, when the wind speed has reached 480 km s~^, at a ra- 
dius of 2.2 X 10^^ pc. Since the cooling time is proportional 
to the third power of the shock velocity, and inversely propor- 
tional to the density, this implies that the cooling is effectively 
30 times higher (the effect of the density is largely cancelled 
due to the fact that the transition happens at a larger radius). 

With the wind velocities mentioned above, we are reach- 
ing the regime where other elements such as Fe, can contribute 
substantially to the cooling. This could mean that the actual 
transition velocity for the solar abundances case is somewhat 
higher. However, our value is in line with what was found 
before, even with cooling curves including these heavier ele- 
ments. Apparently, our omission is not so serious for this case. 
For the [WC] case, where the abundances of the heavier ele- 
ments are mostly 'normal', their contribution to the cooling is 
negligible compared to highly overabundant carbon. 

The expansion velocity of the PN shell is found to be 50% 
higher in the case of the [WC] star, e.g. 21 km s~^ versus 
17 km s^^ att = 2000 years. This is mostly due to the higher 
mass loss rate in the [WC] stellar wind. 

We also ran simulations with a much steeper rise of the fast 
wind velocity, using a = 2.5 in E q. p[ A s was found in previ- 
ous work (ECahn & Breitschwerdt 1990; pwarkadas & Balick 



1998), this does not affect the transition velocities. Obviously, 



the transition times do change. Since the behaviour of the post- 
AGB wind is not well known, one cannot say much about the 
actual transition time, only about the transition velocity. 

This shows that the [WC] character does change the evo- 
lution of the SWB during the early phases of PN formation: 
higher mass loss rates will lead to higher expansion velocities 
and the high C abundance will increase the wind velocity at 
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Table 1 - Abundances used in the various simulations 



Element 


Solar 




[WC] 




WNl 




WNe 




WCl 




WCe 




H 


9.21 X 10~ 


]^ 


9.98 X 10" 


3 


1.51 X 10" 




1.01 X 10" 


-6 


9.89 X 10" 


5 


1.01 X 10" 


4 


He 


7.83 X 10-" 


-2 


7.09 X 10" 


1 


8.38 X 10" 


1 


9.89 X 10" 


-1 


6.43 X 10" 


1 


3.83 X 10" 


1 


C 


3.05 X 10" 


-4 


2.60 X 10" 


1 


1.01 X 10" 


4 


1.01 X 10" 


-4 


2.77 X 10" 


1 


4.64 X 10" 


1 


N 


7.66 X 10" 


-5 


9.98 X 10" 


4 


1.01 X 10" 


2 


1.01 X 10" 


-2 


9.89 X 10" 


5 


1.01 X 10" 


4 





6.22 X 10" 


-4 


2.00 X 10" 


2 


3.03 X 10" 


4 


3.03 X 10" 


-4 


7.91 X 10" 


2 


1.51 X 10" 


1 


Ne 


1.11 X 10" 


^4 


9.98 X 10" 


4 


1.01 X 10" 


4 


1.01 X 10" 


-3 


9.89 X 10" 


4 


1.01 X 10" 


3 



Crowther, private communication 



logm(p/mjj), solar abundances logj^Jp/Tn^^), [WC] abundances 




0.01 O.Oa 0.03 0.04 0.01 O.Oa 0.03 0.04 

Radius <pc) Radius <pc) 



Fig. 2. Colour plot of log^oCp/^H as function of radius and time for the two PN simulations. 



which the SWB changes from momentum- to energy-driven. 
The implications of this and the connection to observations are 
discussed in Sect. 6. 



5. Models for RNe 

For the RNe it is less easy to define one general model. First of 
all, the different types of WR stars have different abundances 
(see Table 1), where one distinguishes between nitrogen-rich 
WN and carbon-rich WC. Largely depending on the effective 
temperature, the subclass '1' for late, i.e. cooler, and 'e' for 
early, i.e. hotter can be added. 

The results from the previous section show that for WC 
abundances, a large effect can be expected. However, for the 
WN stars, the metals are an order of magnitude less abundant, 
helium being by far the most abundant element. Some stars 
are thought to evolve from WN into WC, so the abundances in 
those WR winds will change with time. The other complication 
is the mass loss history of these stars, as pointed out in the 
introduction. 

In order to keep the numerical experiment simple, we run 
simulations of a WR wind running into a Red Super Giant 
(RSG) wind, keeping the abundances in the WR wind con- 
stant. This problem then resembles that of the PNe, except that 
the mass loss rates are higher, and the wind speed is constant 
in time. This way we can study the effects of the H-deficient 
cooling on the bubble evolution. 



We ran five simulations which only differed in the abun- 
dance of the stellar wind, one with solar abundances, and one 
for each type of WR star: WNl, WNe, WCl, WCe. All parame- 
ters are as listed in Table 2. 

Figure ^ shows the evolution of the logarithm of the num- 
ber density in the same way as in Fig. |[ One notices how 
the lowest velocity wind (WNl) model produces a momentum- 
driven RN until t — 550 years, at a radius 7.2 x 10~^ pc, after 
which it switches to energy-driven. The 'solar', WNe, WCe, 
and WCl models are all energy-driven from the start, as is the 
WCe model (not shown). 

Due to the low densities, collisional ionization at the inner 
shock takes some time, and there is transition zone where the 
ionic concentrations are out of equilibrium. This zone is quite 
narrow (2 x 10^^ cm), and its existence does not have a major 
effect on the bubble dynamics or expected X-ray spectrum. 



5. 1. Momentum and energy ratios 



Lamers & Cassinelli (1999) briefly discuss the dynamics of 
SWBs in Chap. 12. There they point out that one can deter- 
mine whether a SWB is momentum-driven or energy-driven by 
estimating the momentum and kinetic energy of the swept-up 
shell, and comparing this to the momentum and kinetic energy 
supplied by the stellar wind during the lifetime of the SWB. 
The ratio of shell momentum to wind momentum is called tt, 
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Table 2 - Input parameters for the RNe models 



Solar WNl WNe WCe WCl 



Msw (Myr-1) 


8 X 10"^ 


8 X 10"^ 


8 X 10"^ 


8 X 10"^ 


8 X IQ- 


Hsw (km s ) 


30 


30 


30 


30 


30 


A^sw (mach number) 


1.1 


1.1 


1.1 


1.1 


1.1 


Mtw (Myr-^) 


1.5 X 10"^ 


3 X 10~^ 


3 X 10'^ 


1.5 X 10"^ 


1 X 10" 


iifw (km s"^) 


2300 


840 


2060 


2300 


1400 


Alfw (mach number) 


20 


20 


20 


20 


20 


ro (cm) 


5 X 10^5 


5 X 10^^ 


5 X 10^^ 


5 X 10^^ 


5 X 10^' 


grid size 


3200 


3200 


3200 


3200 


3200 



logjo(/)/ma), Solar simulation 



1o£io(p/'^h)' ^Cl simulation 



3500 



sooo 



1500 



E- 1000 



500 




15.5B 



4000 



B 3000 



2000 



1000 



0.66 



0.4 



0.1 O.S 0.3 
Radius (pc) 
logj^o(/3/m2), WNe simulation 




15.85 



-0.69 



0.4 



0.1 O.S 0.3 

Radius (pc) 
logj^Q(/)/mg), WNl simulation 




15.76 4000 



0.1 



O.S 0.3 
Radius (pc) 



0.4 




15.97 



0.35 



-0.17 



0.1 



O.S 0.3 
Radius (pc) 



0.4 



Fig. 3. Colour plot of log]^Q(p/TOH) as function of radius and time for four of the five RN simulations. 



the ratio of the kinetic energies is e. 

TT — ; 

M 



(3) 



(4) 



Garcia-Segura & Mac Low (1995 ) derived an analytical de- 
scription for SWBs around WR stars, including expressions for 



the ratios vr and e. However, they assumed that can be taken 
to be zero. While it is true that for WR stars Usw ^ Wfw, often 
^fw^^fw ~ AfswWsw, meaning that one should take the origi- 
nal momentum of the slow wind into account. Also, v^^ can 
be a sizeable fraction of Ushcib which means that the original 
velocity of the slow wind can be a non-negligible fraction of 
the shell velocity. 
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Using a thin-shell approximation, the expansion velocity of 
the shell can be found by considering momentum- and energy 
conservation. For an energy-driven bubble, momentum conser- 
vation means that 



(5) 



where P is the pressure of the hot bubble. This pressure is re- 
lated to input wind energy through energy conservation 



d ( Airr^ 
di 



P \ 1 ■ , o dr 



(6) 



where the radius of the inner shock is assumed to be much 
smaller than the radius of the shell. The mass of the swept up 
shell (Msheii) at a time t, when the shell has reached a radius 



^cxp 



t, is given by 



M 



shell — 



Ms. 



■{r - Vswt) 



(7) 



where Mgw = 47rrQ/9o''^sw, with vq the radius at which the slow 
wind density is pq. 

Equations can be combined into one equation of mo- 
tion for the shell 



4usw dt 



d^ 



r" — {r- Vswty 



dt^ 



1 



(8) 



The solution of this equation is a shell expanding with constant 
velocity. 



^exp,e " 




Vsw-f 












'M 



- 1- Wl-4 




where 




1/3 



2 



(10) 



the solution for inifinitely small v^w, which can be found 
in for example Garcia-Segura & Mac 
Neglecting terms of order (wsw/3w* ^'^ 
mated as 



; Low! ( 



( |1995| ) (Eq. 30). 
can be approxi- 



(11) 



For the momentum conserving case it is easier to consider 
the balance between input momentum and shell momentum. 
The shell now consists of compressed material from both the 
slow and fast winds. The mass of swept up slow wind material 
which we will call A/sh,sw is still given by Eq. ^ the mass of 



the fast wind material compressed into a shell at a position r at 
a time t is given by 



sh.fw 



Vfw 



(12) 



i.e. the amount of mass supplied by the fast wind in a time t, 
minus the material which has not joined the shell yet. 

Conservation of momentum says that the current momen- 
tum of the shell should equal the momentum added by the fast 
and slow winds: 

(A^sh,fw + Msh,lw)Vsh = Af"sh,fwWfw + Msh,{vjVsw ■ (13) 

Substituting Eqs. ^ and |l2| into this expression and some alge- 
bra, shows that the solution is a shell with constant expansion 
velocity 



1'* 



V, 



exp,in — , 

where 



'M, 



fw 



Ms, 



(14) 



(15) 



the solution for Mf^/Ms^ — > (disregarding the mass contri- 
bution of the fast wind), and infinitely small Vsw (disregarding 
the slow wind's momentum input). This was the expression ef- 
fectively used by Garcia-Segura & Mac Low ( 1995| ). 

Only dropping the terms linear or higher in Aft„/Msw and 
Vsw/vfvj, one obtains the approximate solution 

^cxPjin — ^sw ^+,in i (16) 



(see e.g. |Kahn & Bi'eitschwerdt||199(]| , Eq. 16) 

Using these new expressions for the expansion velocities, 
Eqs. ^ and |l^ and assuming a radiative outer shock so that 
i^shcii = Woxp, one can easily derive general expressions for the 
ratios tt and e. Since the expression for the expansion velocity 
in the energy-driven case is quite complicated, we will use the 
approximate expression Eq. ll] to derive a first order approx- 
imation for TT and e in case of a non-zero slow wind velocity. 
This solution is valid as long as AffwWf^ ^ Af. 
obtain 



sw^sw- We thus 



3 



3 (w* 



3 



Y 



For the momentum-driven case, the exact ratios are 



V*.-mVf-w + m 



(17) 
(18) 

(19) 
-(20) 



Taking Vsw = 0, and dropping terms of order w*/wt„, one 
can recover the expressions from parcia-Segura & Mac Lov> 
(|I99^, their Eqs. 83 and 87. 
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Table 3 - Momentum and kinetic energy ratios for RNe models 





Sn1nr 
ooiai 


VV IN 1 


WNp 
VV Lye 


WCp 
VV y^c 


VV 1^1 




2.4 


1.4 


1.5 


1.02 


0.9 




4.1 


2.5 


3.1 


4.1 


4.2 




1.2 


1.1 


1.1 


1.2 


1.3 


e 


0.21 


0.20 


0.16 


0.056 


0.06 




0.43 


0.49 


0.42 


0.43 


0.50 


Em 


0.07 


0.15 


0.09 


0.07 


0.09 



According to Lamers & Cassinelli (1999), RNe are mostly 



momentum-driven. This is mostly based on the work of 



Ireffers & Chu (1982) and Chu (1982), who derived momenta 



and kinetic energies for a set of RNe, and found small values 
for TT and e, typi cally tt '-^ .5 and e ~ 0.01. More recent work. 



for example of |Cappa et al.| (y_996|), gives similar results for 
other nebulae. 



Van Buren (1986) recalculated tt and e for the five RNe 
from Chu (1982). He used the same data, but applied a cor- 
rection for an unseen neutral component. He then finds larger 
values for both tt and e: vr ranges from 1 .2 to 6.8; e ranges from 
0.013 to 0.15. These values would make the SWBs energy- 
conserving rather than momentum-conserving. 

Observationally, one can think of a number of other com- 
phcations when determining vr and e, such as determination of 
the mass and velocity of a clumpy flow, non-sphericity effects, 
unknown mass loss history, estimating the age of the RN, etc. 
As shown above, if the RN is expanding into material from a 
previous slower wind, one would also need an estimate of its 
velocity Wsw 

This raises the question whether observationally derived 
values for tt and e can be used to derive the character of the 
RNe. 

In order to better evaluate the observationally derived val- 
ues, we extract values of tt and e from our simulations. In 
order to check that our procedure works, we determined tt 
and e for two simulations in which the radiative source terms 
were switched off, but which used adiabatic indices 7 = 5/3 
(energy-conserving) and 7 = 1.001 (momentum-conserving), 
respectively. The results from these simulations reproduce the 
analytical values from Eqs. |l7|-^ within 5%. 

In Table 3 we list the values found from the RN simula- 
tions, and the one expected from Eqs. p^7[-|20|. What one sees 
is that even though the SWBs are all energy-driven, in the 
sense that they contain hot shocked fast wind material, the val- 
ues are lower than expected. This means that the partially ra- 
diative character of the SWB reduces the values for tt and e, 
and approach the values for momentum-driven bubbles. For the 
two WC simulations, the value actually drops below the ones 
expected for the momentum-driven case. This is because the 
swept-up material cools below its pre-shock temperature. The 
cooling is more efficient here, due to mixing in of wind material 
into the shell. 

These results shows that the ratios tt and e are of very lim- 
ited use when establishing the character of SWBs: an energy- 
driven bubble may disguise itself as a momentum-driven bub- 
ble. 



6. Discussion 

6.1. PNe 

The effects of the [WC] wind in PNe may seem limited, only 
postponing the transition from momentum- to energy-driven 
bubbles to a later time. However, this later transition may have 
some serious consequences. First of all, during the momentum- 
driven phase the (proto-) PN suffers from the non-li near thin 
shell instability, (see e.g. Dwarkadas & Balick 1998). A later 
transition gives this instability a longer time to work on the 
structure of the nebula. The result could then be a more dis- 
turbed or fragmented nebula. There are some ob servational in - 
dications for more substructure in [WR]-PNe (Gorny 2001). 
More quantitatively, analysis of the emission lines from a num- 
ber of [WR]-PNe shows that these lines are quite wide, some 



times wide enough to hide any splitting of the lines ( Acker et al 



2002). To reproduce these line shapes, the authors invoke a tur- 
bulent velocity component, something which is not needed for 
normal PNe. This could be an effect of the later transition. 

Because of the width of the lines, finding the expansion 
velocity of the nebulae becomes somewhat difficult. Earlier 
results indi cated higher than average expansion velocity for 
[WR]-PNe ( |G6rny|^00l"| ), but [\cker et alj ( |2002| ) dispute this. 
Our models do predict a higher average expansion velocity, 
mostly due to the higher mass loss rates. Of course there will 
be large variations, depending on the precise mass loss history 
of the star, the age of the nebula, and the lines used to measure 
the expansion velocity. 

If the transition from momentum- to energy-driven takes 
longer for [WR]-PNe, there is another interesting effect. 
During the momentum-driven phase, asphericities in the fast 
wind are easily imprinted on the swept-up shell. As the hot bub- 
ble develops, it forms more and more of a buffer between an- 
gular variations in the fast wind and the swept-up shell, and the 
aspherical fast wind has less and less of an effect on the shape 
of the nebula. If PNe acquire most of their asphericity due to 
an aspherical fast wind, one would therefore expect [WR]-PNe, 
which develop a hot bubble later, to be more aspherical. Since 
there seems to be no systematic morphological differences be- 
tween [WR]-PNe and normal PNe ( |G6rny 2001), the conclu- 
sion would be that aspherical post-AGB winds do not play a 
major role in determining the asphericity of PNe. This rather 
qualitative statement will be studied in some more detail in a 
subsequent paper. 

The fact that the [WR]-PNe form a subgroup with clearly 
different mass loss properties may allow further tests of models 
for the formation of PNe in general. For instance, the efficiency 
of the magnetic shaping mechanism of Chevalier & Luo (1994 ) 
is given by 



(21) 



where Bgtar is the magnetic field strength at the surface of 
the star, i?star is the stellar radius, Vrot is the stellar rota- 
tion velocity. If typical field strengths and rotation speeds for 
[WR] central stars are the same as for normal central stars, the 
higher mass loss rates for [WR] stars implies lower values of 
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a, i.e. lesser efficiency of the magnetic shaping. This would 
then give on average rounder [WR]-PNe. The magnitude of the 
effect is about a factor 10. The uncertainties in the other key 
parameters (-Bstar- ^lot), may make it hard to find this effect 
in the observations. However, in all models the formation of a 
PN is basically driven by the post-AGB wind, and that there 
are clear differences between the winds of [WC] and normal 
central stars, a closer comparison of the shapes of [WC]-PNe 
to normal PNe is warranted. 

6.2. RNe 

The numerical models presented in the previous section show 
that it is unlikely that RNe around WR stars are momentum- 
driven. The amount of energy produced in shocking the fast 
wind cannot efficiently be radiated away, not even by the ex- 
tremely metal-rich gas of the WR wind. Unless there are other 
physical mechanisms at work, the RNe cannot be momentum- 
driven. 

A number of physical mechanisms may alter the results 
of our simplified models. Photo-ionization will alter the ion- 
ization state of parts of the RNe. However, the spectrum of 
WR stars is not hard enough to produce the ions which are 
collisionally generated in the hot bubble, and hence it will 
not change our conclusion about the energy-driven character 
of RNe. dumpiness and time variability of the stellar wind 
will produce areas of different temperatures when shocked. 
However, the lower density portions of the wind will cool 
slowly, and are likely to take up the largest volume, so the bub- 
ble will still be filled with a substantial amount of hot gas, mak- 
ing it energy-driven, dumpiness in the environment is likely 
to be more effective in modifying the character of the RNe. If 
these clumps are not effectively accelerated by the outer shock, 
they may end up inside the SWB, and start to evaporate, raising 
the density, and lowering the temperature of the hot shocked 
wind material, a process referred to as 'mass loading'. This 
could conceivably lead to more efficient cooling and a collapse 
of the hot bubble against the inner edge of the shell. Thermal 
conduction, if not inhibited by the presence of a magnetic field, 
can produce a similar result. However, studies of mass loading 



(Pittard et al. 2001 ) show no indication for such a collapse. 



The observed X-ray spectra for RNe show the presence of 
hot gas (Wrigge et al. 1998; Wrigge 1999), although not as 
hot as predicted by the models. This is a strong observational 
indication that RNe are filled with hot gas, and are therefore 
energy-driven. 

The remaining question is how to interpret the observation- 
ally derived values for tt and e. Our simulations show that these 
are of limited value. Even in the idealized case studied here, tt 
and e do not follow the simple division between momentum- 
and energy-driven. The observational problems are also many, 
as was already outlined by Van Buren (1986). In addition the 
complex and fast evolution of these stars and their winds, inval- 
idates the assumption that one can simply integrate over time. 
The momentum and kinetic energy of the shells are an inte- 
gration over the momenta and kinetic energy of the swept-up 
material, which could come from the interstellar medium and 



various previous mass loss phases. Dividing this by a value de- 
rived from assuming the current wind to have been constant in 
time, is unlikely to give sensible values. 

7. Conclusions 

We simulated the effects H-deficient winds have on their 
SWBs, studying simplified cases for PNe and RNe. We find 
that the extreme abundances in the winds of [WC] stars can 
keep their PNe momentum-driven for a longer time. We spec- 
ulate that this leads to more turbulent nebulae and would also 
produce more aspherical PNe if their shape was mostly due to 
an aspherical post-AGB wind. 

For the RNe around massive WR stars, we showed that they 
cannot be momentum-driven, despite earlier claims to the op- 
posite. We pointed out some of the difficulties related to the 
deriving the momentum- versus energy-driven character of the 
RNe using the comparisons of RN momentum and kinetic en- 
ergy to the assumed total input of momentum and kinetic en- 
ergy by the stellar wind. 

The models in this article illustrate the physical effects of 
abundances on the structure of SWBs. To be realistic models 
for PNe and RNe, they should be improved in several ways. As 
shown for instance in Mellema (1994), ionization fronts can 



play an important role in the early evolution of PNe. For the 
RNe adding photo-ionization to the models could help in es- 
timating the amount of neutral material, which is one of the 
unknowns in the tt-e method. Letting their winds evo l ve ac- 
cording to a m ore realistic recipe such as in Melleme ( 1994 ) 
for PNe, and in Garcia-Segura et al. ( 1996a ) for RNe, would be 
another step in the direction of more realistic models. However, 
given that the results in this paper show that for RNe the wind 
abundances p lay on ly a minor role, the mo d els of Garcia- 



Segura et al. ( 1996b ) and 3arcia-Segura et al. ( 1996a ) remain 
largely valid, and only for the case of the [WR]-PNe it makes 
sense to pursue more realistic models. 
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